FROM PARTICLE TO KINETIC AND HYDRODYNAMIC 
DESCRIPTIONS OF FLOCKING 



SEUNG-YEAL HA AND EITAN TADMOR 



Abstract. We discuss the Cucker-Smale's (C-S) particle model for flocking, deriving 
precise conditions for flocking to occur when pairwise interactions are sufficiently strong 
long range. We then derive a Vlasov-type kinetic model for the C-S particle model and 
prove it exhibits time-asymptotic flocking behavior for arbitrary compactly supported 
initial data. Finally, we introduce a hydrodynamic description of flocking based on the C-S 
Vlasov-type kinetic model and prove flocking behavior without closure of higher moments. 
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Collective self-driven motion of self-propelled particles such as flocking of birds and mobile 
agents, schooling of fishes, swarming of bacteria, appears in many context, e.g., biological 
organism [21 [3 El El [13 [HI EH [Ml ES], mobile network [U [H QZB 03] appears in many 
contexts of biological system, mobile and human network [6]. The flocking dynamics 
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of self-propelled particles is important to understand the nature of the aforementioned 
self-propelled particles. The terminology "flocking" represents the phenomenon in which 
self-propelled individuals using only limited environmental information and simple rules, 
organize into an ordered motion (see [23], and it was a subject of biologists [21 [18]. The 
study of flocking mechanism based on mathematical models was first started from the work 
of Viscek et al [25], and was further motivated by the hydrodynamic approach |23| . 

Our starting point is a particle description, proposed recently by Cucker-Smale [3(6], as 
a new simple dynamical system to explain the emergency of flocking mechanism with birds, 
with language evolution in primitive societies etc. The Cucker-Smale's system is different 
from previous flocking models, e.g., [25], in the sense that the collisionless momentum trans- 
fer between particles, {(xi(t) , Vi(t))}fL 1 , is done through a long-range bi-particle interaction 
potential, r(x, y) = r(\x — y\) depending on the distance \x — y\, 



The Cucker-Smale's flocking system (in short C-S system) is reviewed in Section [21 Here 
we revisit the formation of flocking in C-S dynamics in terms of the fluctuations relative to 
the center of mass x c (t) := 1/N^2xi(t). The dynamics of fluctuations makes transparent 
the flocking dynamics. Our main result, summarized in theorem 12. 11 improves [5] for slowly 
decaying interaction potential, r(\x — y\) ~ |x — y\~ 2 @,2(3 < 1. It is shown that flocking 
emerges in the sense that the following two main features occur: (i) the diameter max \x{(t) — 
Xj(t)\ remains uniformly bounded thus defining the "flock", and (ii) the "flock" is traveling 
with a bulk mean velocity which is asymptotically particle- independent, Vi(t) ~ v c := 



When the number of particles is sufficiently large, it is not economical to keep track of 
the motion of each particle through the Cucker-Smale's system. Instead, one is forced to 
study the mean field limit of C-S system and we introduce a kinetic description for flocking, 
in analogy with the Vlasov equation in plasma and astrophysics. In Section [3J we present 
a Vlasov type mean field model, which is derived from the C-S system using the BBGKY 
hierarchy in statistical mechanics. The formal derivation, carried in Section [31 follows by 
taking the limit of an iV-particle interacting system consisting of self-propelled particles 
governed by C-S flocking dynamics. To this end, let / = f(x,v,t) denote the one-particle 
distribution function of such particles positioned at (x, t) G M rf x M + with a velocity vel^: 
the dynamics of the distribution function / is determined by 



dictated by a prescribed interparticle interaction kernel, r = r(x, y). We refer to Degond and 
Motsch [71 [HI [9] for recent kinetic description of Vicsek type model of flocking. The dynamics 
of particle trajectories of the proposed kinetic description of flocking is analyzed in Section 
13.21 in Section 13.31 we prove the global existence of smooth solutions to the kinetic model 
with arbitrary smooth compactly supported initial data. In Section [H we show that the 
kinetic model reveals the time-asymptotic flocking behavior when the bounded interparticle 
interaction rate has a sufficiently strong long range. Our results are summarized in the 



Tt Vi ^=N ^ r(x i (t),x :j (t))(v 1 (t)-v l (t)). 



l/ivxxo). 





r{x,y){v* - v)f{x,v,t)f(y,v*,t)dv*dy, 
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main theorem l4.lt proving the decay of energy fluctuations, A [/](£), around the mean bulk 
velocity, u c , 

r , , s f ,o , v , . \-o2dV fix, v,{)dvdx , , 

A / (t) := / t>-Uc a /foM)dt;cfa;, tt c t = f J, \L 3 = ^ » ■ 
J R 2d J R2d f(x,v,t)dvdx 

Flocking is proved for the restricted range, 2/3 < 1/2, realized by the asymptotic decay 
estimate, A[f](t) — ► as i — > oo. 

In Section [5] we turn our attention to the hydrodynamic description of flocking, furnished 
by moments of the kinetic distribution function. We study the dynamics of the resulting 
system of balanced laws related to the moments of Vlasov model. Despite the lack of closure, 
we present a fundamental estimate which enables to conclude the flocking mechanism at the 
macroscopic hydrodynamic scales. Theorem 15.11 states that the energy-related functional, 

r(t) 

r (*) : = / ( n\ u ( x ) -u(y)\ 2 + e(x) + e(y))p(x)p(y)dydx. 

J R 2d \2 J 

decays provided the interparticle interaction, <p(s) = inf r(x(s), y(s)) decays slowly enough 
so that its primitive, $(t), diverges. This in turn in related to the increase of entropy 



flog(f)dxdv>0, 

at J R 2d 



as particles with increasingly highly correlated velocities flock towards particle-independent 
bulk velocity. 

2. A PARTICLE DESCRIPTION OF FLOCKING 

2.1. The Cucker-Smale model. In this section, we briefly review the Cucker-Smale's 
flocking system in [5, 6, 22j, which manifests the time-asymptotic flocking behavior of many 
particle systems. We reinterpret the C-S system in terms of fluctuations relative to the 
center of mass coordinates, which enables us to simplify and sharpen the derivation of 
sufficient conditions for flocking to occur. 

Consider an A^-particle interacting system consisting of identical particles with mass m to 
be assumed to be unity. Let [xi(t) , Vi(t)] £ M. 2d be the phase space position of an i-particle. 
The Cucker-Smale dynamical system [6] takes the form 

(2-1) j t xi(t) = Vi (t), i^(i) = A rixiitlxjWivj^-Viit)). 

Here A is a positive constant, and r(x, y) is a symmetric, bi-particle interaction kernel, 
(2.2) r(x,y) = r(y,x) < A. 

To discuss the time asymptotic flocking behavior, we will restrict our attention to inter- 
particle interactions which are decreasing functions of the distance^, 

(2.3a) r(x, y) = r(\x — y\), r(-) isdecreasing. 

A prototype example is the interaction kernel with a polynomial decay of order 2/3, [5], 

(2.3b) r(\x-y\)> A > 0. 

(1 + \x — y\ )' 



^To emphasize this point, we therefore continue to refer to general symmetric kernels, r(x,y) = r(y,x) 
whenever translation invariance is not necessary. 
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We note in passing that only a lower-bound of the interaction kernel matters. 

For notational simplicity we often omit i-dependence from the particle identification, 
abbreviating x^ = Xi(t) and Vj = Vi(t). 

Let rrij(t), j = 0, 1, 2 denote the moments. 

N N N 

m :=^l = iV, m x (t) :=^2vi(t), m 2 (t) := ^ \vi(t)\ 2 . 

i=l i=l i=l 

Regarding the dynamics of these moments, we have the following estimate. 

Proposition 2.1. Let (xi(t),Vi(t)) be the solution to the C-S system (|2.ip . (|2,2p . Then the 
following estimates hold. 

(2.4a) -mi(i) = 0. 

(2.4b) ^t 1712 ^ = ~Jj Y r(x h x j )\v j -Vi\ 2 . 

(2.4c) m 2 (t) > m 2 (0)e- 2XAt + l^lM! ( 1 _ 6 "2AaA 

mo V / 

Remark 2.1. Proposition \2.1\ tells us that although the kinetic energy mzit) is mono- 
tonically decreasing, \2.^b ), it has the following nonzero lower bound if initial momentum 
mi(0)^0, 

(2.5) m 2 (t) > l^lM. 

m 

Flocking occurs when equality takes place in the Cauchy-Schwarz inequality Ii2.5\) . 

Proof. Conservation of momentum in (|2.4ap follows from the symmetry r(xi,Xj) = r(xj,Xi), 
for 

d/A \ A 



Y Vi ) = N Y r ( X i, X j)( V 3 ~ V i) = °- 



dt\^ J N 

i=l l<i,j<N 

Moreover, symmetry also implies 

^ r(xi, Xj)vi ■ (vi - vj) = - ^2 r(xi,Xj)vj ■ («» - Vj) 
l<i,j<N l<i,j<N 

and hence the energy dissipation (|2.4bp follows 

J t (^Z\ v i\ 2 ) = ~Jf Y rix^x^Vi-ivi-Vj) 
i=l l<i,j<N 

= Y r ( x i, x j) v j ■ ( v i - v j) = ~jj Y r(x h Xj)\vi-Vj\ 2 . 

l<i,j<N l<i,j<N 
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Finally, to prove (|2.4c|) . we use the energy dissipation in (|2.4b|) . the fact that r(xi,Xj) < A 
and the conservation of momentum in (|2.4aP to find 

_ m2 (i) = _A ^ r(xi,xj)\vi -Vj\ 2 > -— A Y \vi-Vj\ 2 

l<i>i<iV l<i,j<N 

= -2XA( m2 (t) - = -2XA( m2 (t) |mi( ° )|2 



N ) V mo 

Gronwall's lemma yields (|2.4cp . □ 

2.2. Asymptotic behavior of fluctuations — flocking. We now turn to study the 
asymptotic time behavior of solutions to C-S system (|2,ip . (|2.3ap . To this end, we introduce 
a center of mass system (x c (t),v c (t)), 

1 N 1 N 

i=l i=l 

Then, thanks to conservation of momentum, the velocity v c is constant in t, and the tra- 
jectory of center of mass x c is a straight line: 

v c {t) = u c (0), x c (t) = x c (0) + tv c (0). 

Observe that the fluctuations around the center of mass, 

Xi(t) ' ^ Xi{t) - x c (t), Vi(t) i ^ - u c (i), 

satisfy the same C-S system (|2,ip . (|2,3ap : it is here that we take into account the fact that 
the interparticle kernel depends on the distance, r(x,y) = r(\x — y\). We shall show that 
under appropriate conditions, flocking occurs in the sense that these fluctuations decay in 
time. Thus, the time-asymptotic dynamics of C-S solutions emerges as a linear movement 
with a fixed velocity dictated by the coordinates of center of mass. 

To proceed, we introduce the two auxiliary functions which measure the fluctuations of 
the fluctuations around their center of mass, 

X(t):= £ \x t {t)-x c {t)\ 2 , V(t):= £ \ Vi {t) - v c (0)\ 2 , 
l<i<N l<i,j<N 

subject to initial conditions (Xo,Vq) = (X(0),V(0)). The flocking behavior will depend in 
an essential way on the behavior of the minimal value of the interparticle interaction at 
time t, 

(2.6) ip(t) := mm r(xi(t),Xj(t)). 

1<i,j<N 

We begin with the fluctuations of velocities. 

Lemma 2.1. [Fluctuations of velocities]. Let (xi(t),Vi(t)) be the solution of the system 
()2.1|) . (l2.3ap . Then we have 

V(t) < V e~ 2xm , $(t) := I cp(r)dT. 

Jo 

Proof. We invoke (|2.4bp with Vi(t) — v c (t) replacing Vi{t) to find that 

d N A 
( 2 -7) ~ Vc \ 2 = ~N Yl r{x i ,x j )\v j -Vi\ 2 . 

i=l l<i,j<N 
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Since ~~ v c) = 0, we have Yli<ij<N \ v i ~ v j\ = 2NV(t), and the result follows from 

Gronwall's integration of 

Jt V{t) = ~Jl ^ rixux^vj-vil 2 <-2\<p(t)V(t). 

□ 

Remark 2.2. Lemma \2. 1\ implies the sufficient condition for flocking is that the interparticle 
interaction potential decays sufficiently slow, so that its primitive, 3>(i), diverges, i.e., 

(2.8) if lim $(t) = / (p{r)dT = oo then lim \vi(t) - v c \ = 0, i = l,---,N. 

t^OO Jq t^OO 

The answer whether <p(t) decays sufficiently slow to enforce flocking depends on the 
variance of positions Xi(t). 

Lemma 2.2. [Fluctuations of positions]. Let (xi(t),Vi(t)) be the solution of the system 
(|2.ip . (|2.3ap . Then we have 

t 2 

X(t)<2X + V -, t>0. 
Proof. We use Cauchy-Schwartz's inequality to see 

N N 



d_ 

dt 



£ 

i=l 



Xi Xr 



dX'i dx £ 
~dt dt 



i=l ^ 
N 

}Xxj - %c) ■ (Vi - v c ) < 



8=1 

d 



N 



Xf X r- 



N 



Vi - V r 



Using Lemma Owe obtain, ^X(t) < y/V(t)y/X(t) < y/%e~ x<s>(f) y/X(t), and the solu- 
tion of this differential inequality yields 



(2.9) 



Vn 

X(t) < 2X + y 



< 2X + Vo-. 



□ 



As a corollary of Lemma 12.21 we now obtain the desired lower bound for (p(t). 
Corollary 2.1. Let (xi(t) , Vi(t)) be the solutions to (|2.ip . (|2.3ap . Then cp(t) satisfies 



<p(t) > r(y/2X(t)) > r(V4X + y t 2 ). 
Proof. We use Lemma [2721 to see that for each i, j £ {1, ■ ■ ■ , N}, 

\ Xi - xj\ 2 < 2(\ Xi - x c \ 2 + \xj - x c \ 2 ) < 2X(t) < 4X + V t 2 . 



Since r(-) is decreasing, we have ip(t) = min r(xi(t),Xj(t)) > r (a/ AXq + Vot 2 ) . 

1 <i.j< X 



□ 



The asymptotic flocking now depends on the specific decay of the interparticle interaction 
r(-). As an example, consider the C-S system with the 2(3 interaction (|2.ip . (|2.3bp . where 



(2.10) (p(t) > A(l + 4X + V t 2 y p > Am(l + t 



-2/3 



Ki := ( max{l + 4X0,^0} 



-P 
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We conclude that the divergence of $(t) = f (p(r)dr and hence, by (12. 8p that flocking 
occurs, for 2(3 < 1. This recovers the Cucker-Smale result [U E|. Below we improve the 
Cucker-Smale result proving unconditional flocking result for (3 = 1/2. 

Theorem 2.1. Let (xi(t) , Vi(t)) be the solutions to (J2~7lj) . (|2.3bl) with V > 0. Then the 
following holds. 

(i) There exist a positive constant, C 2 (depending only on k,\,A and (3 as specified in 
(|2.12p below), such that 

(2.11a) \xi{t) - x c \ < \ Xi (0) - x c \ + C 2 . 

(ii) There exists constants, Ki > 0, i = 1,2, such that 

e -XA K2 t^ /3e[0,i), k 2 := (1 + 4X + SCf)-/ 3 , 



(l + t)~ XAK \ (3 = ± Kl = ( max{l + 4X ,y } N 



(2.11b) \ Vi (t)-v c \ < y%x { 



Remark 2.3. Theorem \2.1\ shows the two main features of flocking occur with the 2(3- 
inter action potential, 2(3 < 1, namely, the diameter max \xi{t) — Xj(t)\ remains uniformly 
bounded thus defining the traveling "flock" with velocity which is asymptotically particle- 
independent, Vi{t) v c . 

Proof. We begin with the case < (3 < ^. To get the optimal exponential convergence rate, 
we employ a bootstrapping argument in three steps. 

Step 1. We first obtain a weak integrable decay rate for \vi — v c \. Using (|2.10p we find 
(2.12a) - j\(T)dr <-A Kl ^(l + rr^dT^C^l-il+t) 1 - 2 ^, < (3 < X -. 

The above estimate together with Lemma 12. II yield V(t) < Voe™ 2 ^ 1 "^ 1 2/3 . 

Step 2. Next, we improve lemma l2T2| observing that for (3 < 1/2, the position X(t) remains 
uniformly bounded in time. Indeed, we have, 

Xi(t) = Xi(0) + / Vi(r)dT, x c (t) = x c (0) + / v c dr. 

Jo Jo 

Time integr ability of \vi(t) — v c (t)\ then yields (|2.11ap . 

\xi(t) -x c (t)\ < \xi(0) -x c (0)\ + [ \vi{r)-v c \dr 

Jo 

roo 

< \ Xi (0) - x e (0)\ + / e- c ^ l+t ^- 2f3 dt < |xi(0) - x c (0)| + C7 2 , 
Jo 

where 

(2.12b) C 2 < / e-°^ 1+t ^dt < oo, 0<(3<-. 

Jo 2 

Step 3. The uniform bound of X(t) implies an improved estimate for the interparticle 
interaction <p(t): corollary 12.11 implies 



<p(t) > r(^/2X~ Q ) > A(l + 4X + 8C|) -/3 = Ak 2 , 
which in turn, using lemma [2TT| yields the optimal exponential convergence rate (I2.11bl) 
Ht) - v c (t)\ 2 < V(t) < V e~ 2X ^ < e- 2XA ^. 
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It remains to deal with the case j3 = \. Here, we have 

-2A*(t) < -2AA(^max{l + 4X , F }) 2 J (l + r^dr = -2\Ak\ ln(l + i), 

which in turn implies (I2.11bl) . V(t) < V r e" 2A * ( * ) < + ^~ 2Aj4ki . □ 
Remarks. 

(1) Consider the borderline case (3 = \ with initial configuration satisfying 



\k\ > 1, i.e., \J max{l + AXq, Vq} < A; 

then the same bootstrapping argument used for (3 6 [0, h) gives the exponential 
convergence: 

\vi(t)-v c \ < v^e"^ 2 '. 

(2) Flocking occurs even if (5 > i, but only for special initial configurations. Sufficient 
flocking conditions for such initial profiles is presented in [5]. 

3. From particle to kinetic description of flocking 

3.1. Derivation of a mean- field model. We assume that the number of particles involved 
in the C-S model (|2.ip . (|2.2p is large enough that it becomes meaningful to observe the N- 
particle distribution function, 

(3.1a) f N = f N (x 1 ,v 1 ,...,x N ,v N ,t), (x i)Vi ) £R d xR d . 

Since particles are indistinguishable, the probability density / = f N (■) is symmetric in its 
phase-space arguments, 

(3.1b) f (• • • , Xi, Vi, ■ ■ ■ Xj, Vj, ■ ■ ■ ,t) = f (• • • , Xj , Vj , • • • X{, Vi, ■ ■ ■ ,t), 

so we can 'probe' / by any of its N pairs of phase- variables. Let f N (-,-,t) denote the 
marginal distribution 



/ 



J W 2d(N-l) 

The formal derivation of a kinetic description for the C-S particle system ()2.ip . ()2.2p is 
carried out below using the BBGKY hierarchy, e.g., Q3JJ [5D], based on the Liouville 
equation, [15] 



N , N N 

(3.2) d t f N + ^ E V ^ • (»-'---'-^0 " v ^f N ) = 0- 

i=l i=l j=l 

To this end, one study the marginal distribution f N (x\,vi,t) by integration of (13. 2p with 
respect to dx^dv^ = dv2dx2 • • • dvj^dxjy (to simplify the notations, we now suppress the 
time-dependence whenever it is clear by the context, denoting f N {x\,V\, ■ ■ ■ ,XN,VN,t) = 
f N (xi,v\, ■ ■ ■ ,xn,vn)). Since f N (-, •) is rapidly decaying at infinity, the transport term in 
(|3.2p amounts to 

N 

Nj„ t . _ ... w fN, 



(3-3) / X>" V ^ 



dx-dv- = v\ ■ V Xl f (xi,v\). 
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The corresponding integration of the forcing term in (|3.2p . yields 
A N f N 

— V / • (r(xi,Xj)(vj -Vi)f N )dx-dv- 

/ " ( r ( x i' x j)( v j - vi)f N )dx-dv-. 

7 M 2d(jv-i) f-f v / 



A^ 



2<j<N 



But the symmetry of (|3.1b|) . implies that the integrals being summed above are the 
same for j = 2, 3 . . . , N. Consequently, it will suffice to consider j = 2: 



4$Z / J2 V n- (r(x l ,x j )(v j -Vi)f N )dx-dv^ 

(3.4) = ^-(N-l)f r(x 1 ,x 2 )V Vl ■ ((v 2 - v 1 )f N )dx 2 dv 2 - ■ ■ dx N dv N 

JS J R 2d(N-l) V / 

A- — jV^-f / r(xi,x 2 )(v 2 - v 1 )g N dx 2 dv 2 ). 



Here g N is the two-particle marginal function 

JV/ „, „ „, .\ . / /-A 



5 (zi,Ui,x 2 ,U2,*) := / / dx 3 dv 3 - ■ ■ dx N dv N . 



Thus, in view of (|3.3p and (|3.4p . marginal integration of (|3.2p over implies that 

the one-particle density function, f N (x\,vi,t), satisfies 

d t f N + vi • V X J N + V«! • ( jf m rfxi, a? a )(«2 - wi^cteadt*) = 0. 

We now take the mean-field limit N — > oo: we end up with the one- and two-particle limiting 
densities, / := limjv-»oo f N (%i,Vi) and 9 := limjv->oo g N (%i, vi, x 2 , v 2 ), which satisfy 

(3.5) d t f + vi ■ V xi f + XV Vl ■ ( / r(xi,x 2 )(v 2 -v 1 )gdx 2 dv 2 = 0). 

To close the above equation we make the "molecular chaos" assumption about the indepen- 
dence of the two-point particle distribution, 

g(x 1 ,vi,x 2 ,v 2 ,t) = f(x 1 ,v 1 ,t)f(x 2 ,v 2 ,t); 

Relabel, (xi,vi) i— ► (x,v) and (x 2 ,v 2 ) i— > (y,v*). We conclude that the one-particle distri- 
bution function f(x,v,t) satisfies the Vlasov-type mean-field model, 

(3.6a) d t f + v V x f + \V V • Q(f, f) = 0, 

(3.6b) Q(f,f)(x,v,t):= r(x,y)(v*-v)f(x,v,t)f(y,v*,t)dv*dy. 

JR 2d 

Here, Q(f, /) is the quadratic interaction which can be expressed in the equivalent form 
(3.6c) Q(f,f)(x,v,t)=fL[f], L[f](x,v,t) := r(x,y)(v* - v)f(y,v*,t)dv*dy. 
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3.2. A priori estimates. We begin our study with a series of a priori estimates on the 
solution of the mean-field model fl3.6|) . and the growth rate of the x and ^-support of /. We 
first set 

M0-=^ M0~Zi i = l,-..,d, and Vd+i(0 := l£| 2 - 
Let / be a classical solution to (|3.6p with a rapid decay in phase space M. 2d . A straightforward 
integration of (13, 6p yields 

(3.7a) ^- / tpi(v)f(x,v)dvdx = [ V v ipi(v) ■ Q(f, f)dvdx, 

at j R 2d j R 2d 

(3.7b) ± «,)/(,,„)** = / V.(* W - .)/(,,,)**. 



Using (|3.7p we obtain 

Proposition 3.1. Let f be a classical solutions decaying fast enough at infinity in phase 
space. Then following macroscopic quantities associated with f, satisfy 

(3.8a) — / vf(x,v)dxdv = 0; 
at J R 2d 

(3.8b) ±[ W V(«,.)dW,~/' r( x, !/ )|„-„»|V(x,»)/fe,».)*,^; 

dt jR2d J R 4d 



(II 

Proof. Equality (|3.8a|) follows from (|3.7ap with ipi( v ) = v ii 

V v ipi(v) ■ Q(f,f)dvdx = / r(x,y)(v*i — Vi)f(y,v*)f(x,v)dv*dydvdx = 0. 



The last integral vanishes due to antisymmetry of the integrand, realized by the interchange 
of variables (x,v) <-> (y,v*). The statement of (j3.8bj) follows from (|3.7bp with i^d+i{v) = 
\v | 2 , and observing that 



V^d+l (v) ■ Q{f, f)dvdx = 2 1 v ■ Q(f, f)dvdx 

JM. 2d 

= 2 / r(x, y)v ■ (u* — v)f(y, v*)f(x, v)dv*dydvdx 

JK. id 



r(x,y)v* ■ (v* — v)f(y,v*)f(x,v)dv if dydvdx 

id 

r(x, y)\v — v*\ 2 f(x, v)f(y, v*)dv*dydvdx. 
Finally, we note the two identities, f v ~ x v ■ V x f = -v ■ V x f p , and 

f p -^v ■ Q(f, f) = V v - (^y^-) + (l - ~) (v„ • L[f]) r, Q(f, f) = fL[f\. 

Integration of (13. 6p against f p ~ 1 then yields 
d 
dt 



f p dvdx 



-P I P' 1 (v • V x f + V, • Q(/, /)) dvdx 
-(p-1)/ (v,-L[/])/«x 
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= -d(p-l) r(x,y)f(y,v*)f p (x,v)dv !l: dvdydx. 

JM. 4d 

□ 

Let / be a classical kinetic solution of (|3.6p . The statement of (|3.8cp shows that its 
L 1 (dxdv)noim, the total macroscopic mass is conserved in time (while according to (13.8cp . 
higher L p (dxdv)-norms of / decay in time), 

(3.9a) M Q (t) := / f(x,v,t)dxdv = M . 

Jw. 2d 

Similarly, (|3.8ap tells us that the total macroscopic momentum is conserved in time, 
(3.9b) Mi(t):= vf(x,v,t)dxdv = Mt. 



Finally, (|3.8b|) tells us that the total amount of macroscopic energy is non-increasing in 
time, 

(3.9c) M 2 {t):=f \v\ 2 f(x,v,t)dvdx<M 2 (0). 

JE. 2d 

Here, M.q := Mo(0), M\ ■= A^i(O) and ^2(0) denote, respectively, the initial amounts 
of mass, momentum and energy at t = 0. Next, we turn to the following a priori bound on 
the kinetic velocity. 

Lemma 3.1. Let [x(t),v(t)] be the particle trajectory issued from (x,v) £ suppr xv \fo at 
time 0. Then the i-component of velocity trajectory, Vi(i) = Vi(t;0, x,v), i = 1, ••• ,d, 
satisfies 

Vi (t) S ( Vi (0)e- XAMot - A(i - e ~ XAMot ), Vi (0)e- XM ^ + XAJ f e - XM ^~^ds 
^ Mo Jo 

Here, 3>(i) := f Q t ip(s)ds, M.q = ||/o||zi ^ s the initial total mass and J := ^/^AoM2■ 

Proof. For given (x,v) £ swpp xv fo, we set x(s) = x(s;0,x,v) and v(s) = v(s;0,x,v). Note 
that for each i = 1, • • • , d, we have 



(3.10) Li[f(x(t),v(t),t)]= r(x(t),y)(y*i(t) -Vi)f(y,v*)dv*dy 

JR 2d 

= Vi (t)+ r{x(t),y)v*if{y,v*)dv*dy-[l r(x(t),y)f(y,v*)dv Jf dy). 

JB. 2d v Jm 2d ' 

Lower and upper bounds for the kinetic velocities are obtained in terms of the estimates, 
<f(t)M < / r(x(t),y)f(y,v*)dv*dy < AM , 



r(x,y)v*f(y,v*)dv*dy < A^\\fo\\ L i v ^\\\v\ 2 f \\ L i v = AJ . 

It then follows from (|3.10p that 

-XAMoVi(s) - XAJo < ^(t) = \Li[f(x{t),v(t),t)] < -Xtp(t)M Vi(t) + XAJ , 

dt 

and the desired result follows by Gronwall's integration. □ 
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Remark 3.1. Let S2(i) denote the v-projection o/supp/(-, t), 

(3.11) ft(t) := {d e l d : 3 (x, u) G R M sttc/t t/iat /(x, v, i) / 0}. 

Lemma \3. 1\ shows that if fo(x,-) is compactly supported, then supp/(x,-,i) remains finite, 
with a weak growth estimate for the velocity trajectory, 

(3.12) \vi(t)\ < max{r/o + -r^-,Vo + XJot\ < r?0 + 77- + XJ t, 770 := max 

1. Mq J .Mo ugn(o) 

3.3. Global existence of classical solutions. In this section we develop a global exis- 
tence theory for classical solutions of the Vlasov-type flocking equation, 

(3.13a) d t f + v-V x f + XV v -[fL[f])=0, x,v£R d ,t>0, 

f A 

(3.13b) L[f](x,v,t)= r(x,y)(v* - v)f(y,v*,t)dv*dy, r(x,y) = — — 

Jr™ {l + \x- y\ 2 )P 

subject to initial datum 

(3.13c) f(x,v,0)=f (x,v). 

We begin by noting that the kinetic solution / remains uniformly bounded. To this end, 
rewrite the mean-field model f|3. 13|) in a 'non-conservative' form, 

(3.14) dt f + vV x f + XL[f}-V v f = -XfV v -L[f], x,v£R d ,t>0. 

Consider the particle trajectories, [x(t),v(t)] = [x(t; to, xq, vq), v(t; to, xo, vq)], passing through 
(x ,vq) G R d x R d at time t G R + , 

(3.15) j t x(t) = v(t), j t v(t) = XL[f(x(t),v(t),t)]. 

Noting that -V„ • L[f] = d J R2d r(x,y)f(y,v*,t)dv*dy, we find 

||V«.L[/]|U« < dA\\f\\ L i iV = dAM , 
which implies that the following inequality holds along the particle trajectories, 

j t f{x{t),v{t),t) < XdAMof{x{t),v{t),t). 

It follows that as long as initial data fo has a finite mass, there will be no finite time blow-up 
for /(•,*), 

(3-16) \\f{t)\\L~ v <e UAM ^\\fo\\L^. 

Next, we turn to study the smoothness of f(-,t). Since the local existence theory will be 
followed from the standard fixed point argument, e.g., 3], we only obtain a priori C 1 -norm 
bound of / to conclude a global existence of classical solutions. 

Theorem 3.1. Consider the flocking kinetic model $3.13]) . Suppose that the initial datum 
fo G (C 1 n W 1 ' 00 )^ 2 *) satisfies 

(1) Initial datum is compactly supported in the phase space, supp( XiV \fo(-) is bounded, 
and in particular, 17(0) C B vo (0). 

(2) Initial datum is C 1 -regular and bounded: 

£ n v *v£/oik~ < 00. 

0<\a\+\p\<l 
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Then, for any T 6 (0, oo), there exists a unique classical solution f G C7 1 ([0,T) x M. 2d ). 

Proof. We express the non-conservative kinetic model (13. 14[) in terms of the nonlinear trans- 
port operator d t + v ■ V x + XL[f] • V v , 

(3.18) Tf = -XfV v -L[f], T:=d t + vV x + XL[f]-V v . 

We claim that there exist (possibly different) positive constants, C = C(d, X, Mo, Jo) > 0, 
such that 

(3.19a) |T(/)| <C\f\, 

(3.19b) \T(d x J)\ < c(\f\ + ( V (t) + 1)|V„/| + \d x J\), n(t) := max >|, 

(3.19c) \T{d v J)\<c(\d x J\ + \V v f\). 

To verify these inequalities, observe that by (13.181) 

= -XfV v ■ L[f] = Xd I r(x,y)f(y,v*,t)dv*dy 

JR 2d 

and (l3TT9aD follows with C := XdAM > X\\V V ■ L[f](t)\\ L «> v . 
Next, differentiating (|3.18p we obtain 

T(d x J) = -X(d Xt L[f\) ■ V,/ - X(d Xi V v ■ L[f])f - X(V V ■ L[f})d x J. 
Straightforward calculation yields 

dxiL[f] = - [ t^^ X% \2\3+i ( v * -v)f(y,v*,t)dv*dy, 
J R 2d (l + \x- y\ z ) p+1 

and since the variation of the relevant kinetic velocities at time t does not exceed \v — u* | < 
2n(t), we find \\d Xt L[f](t)\\ < A/3Arj(t)M ; similarly, 

d Xi V v -L[f)= [ 2 ^*fo - /( t)M _ \\ d Vv . L [f](t)\\ L < 2 pdAMo. 
J R 2d (1 + \x- y\^) p+1 

We conclude that (|3.19bl) holds with, say, C = XdAM Q (l + 2(3 + 4/3ry(t)). 

Finally, we differentiate (|3.18p with respect to Vi (noting that d Vi V v ■ L[f] = 0) 

T(d v J) = -Xd x J - X(8 Vi L[f}) ■ V v f - X(V V ■ L[f])d v J; 

Straightforward calculation then yields, 

dvMf] = - r(x,y)f(y,v*,t)dv*dy i-> \\d Vi L[f](t)\\ L ^ < M , 

and (l3TT9cl) follows with C = X + X{d + l)AM . 
Now, let J~(t) measure the VF 1,00 -norm of f(-,t) 

m ■■= E n v ^/(*)ii^- 

o<M+|/3|<l 

The inequalities (|3.19|) imply 

j t Ht)<Ht) + i)T{t). 

Lemma [37TI (see remark [3TT]) . tells us that r](t) < rjo + 1, and we end up with the energy 
bound 

F(t) < H0)e C(t+t2) , C = C( m ,Mo,Jo,P,d,A). 
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Equipped with this a priori PI/ 1 ' 00 estimate, standard continuation principle yields a global 
extension of local classical solutions. □ 

Remarks. 

(1) The above a priori estimate need not be optimal. Since we used a rough estimate 
(137L21) for the size of fi(t), 

max |u| < 770 + t, 

vEU(t) 

we end with the quadratic exponential growth, e c '( t +* 2 ) . An optimal bound, however, 
could be e . Of course, one cannot expect a uniform bound for C 1 -norm, because 
the one-particle distribution function may grow exponentially along the particle 
trajectory (see (|3.16p ). 

(2) The global existence of classical solution can be improved for more general kernels. 

(3) For related works on kinetic granular type dissipative systems, we refer to [12 H 13 | H3]. 



4. Time-asymptotic behavior of kinetic flocking 

In this section, we present the time-asymptotic nocking behavior of the kinetic model for 
flocking (I3.13p . As in the case with particle description discussed in Section [2j we will show 
that the velocity of particles will contracted to the mean bulk velocity u c , which corresponds 
to the velocity at the center of mass: 

u c{t) '■= -r7~ I vf(x,v,t)dvdx, u c (t) = u c (0). 



Mo 

We recall that the energy decay in (|3.8bj) 
(4.1) j t M 2 (t)<0. 

We note that unlike granular flows, for example, e.g. [3], the energy decay (|4.ip does not 
drive the energy to zero: if the initial momentum Mi ^ 0, then the kinetic energy M2(t) 
has a nonzero lower bound, in analogy with the discrete case, consult remark \2.1l Indeed, 
since the total mass and momentum are conserved, Mi{t) = Mi, i = 0, 1, (|3.8bp implies 



d_ 

dt 



M 2 (t)>-A / \v-v 3f \ 2 f(x,v)f(y,v*)dv*dydvdx = -2AMoM 2 (t) + 2A\M 1 
Jm 2d 



and Gronwall's lemma yields the following kinetic analog of (|2.4cj) . (j2.5|) 

(4.2) M 2 (t) > M 2 (0)e- 2AMot + 1^(1 - e~ 2AMot ) > 

Mo Mo 

Thus, energy decay by itself does not assert flocking. As in the particle description, the 
emergence of the time-asymptotic flocking behavior depends on the sufficiently slow decay 
rate of the interparticle interaction tp(s) = r(x(s),y(s)). To this end, we let S denote the 
x-projection of supp/(-,i) 

(4.3) E(t) := {x e R d : 3 (x, v) e R 2d such that f{x, v, t) ^ 0}, 
and denote its initial size, Co '■= niax^gSo \ x \- 
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Lemma 4.1. Let f be a global classical solution to ()3.13j) . Then, there exists a k$ > such 
that <p(t) = M( Xj y) e s( t ) r(x, y) satisfies 

(4.4) <p(t) > Ak~ 13 (I + t 2 + t 4 )"^. 
The constant n 3 is given by 

k 3 := max { 1 + 12C 2 , 12 + A) 2 , 3A 2 J 2 } . 

Proof. Let (x,f) G supp^/o- It follows from remark [37T1 that 

(4.5) |«i(t)| <% + 77- + AJ t, 
and hence 

AJ t 2 



\xi(t;0,x,v)\ < \xi\ + J \vi(s;0,x,v)\ds <Co+ (r] + -^-^t + 



2 

This gives an estimate on the size of x-support £(t) of /. For x,y £ S(t), 

C ( t ) < 1 + \ x _ y|2 < 1 + 2 (| X | 2 + |y| 2 ) 



< 1 + 4 



Jo \ , AJ t 2 i2 



^ + ^ + M-J t+ 2 
Jo ^ 2 



< 1 + 12C 2 + 12 % + t^T * 2 + 3A 2 J 2 t 4 

< K 3 (l+t 2 + t 4 ), 

and dO]) follows, <p(t) > r(C(t)). □ 

Let v — u c denote the fluctuation (or peculiar) kinetic velocity. We will quantify the 
emergence of the time-asymptotic flocking behavior in term of the corresponding energy 
fluctuation 

(4.6) A[/(i)] := / \v -u c \ 2 f{x,v,t)dvdx. 

The time-evolution estimate of A[f(t)], will depends on the decay rate of <p{t). 

Let / be a classical solution of f)3. 13H with compact support in x and v. Direct calculation 
implies 

d . r „,.„ r , |2 , 



-A[/(i)] = / |«-u c | d t f(x,v)dvdx 



= — I \v — u c \ 2 v ■ V x fdvdx — A / |u — u c | 2 V„ • f)dvdx =: Zi + X2- 
'R 2d 'r m 

The first term on the right vanishes by the divergence theorem 



T\ 



v — u c \ 2 v ■ V x fdvdx = — / div x ( \v — u c \ 2 vf)dvdx = 0. 
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The second term is simplified as follows. 

1 2 = 2X [ (v-u c )-Q(f,f)dvdx 
Jm. 2d 

= — 2A J r(x,y)(v - u c ) ■ (v - v*) f '(y ,v*) f '(x ,v)dv*dvdydx 
= — 2A / r(x, y)v ■ (v — v*)f(x, v)dv*dvdydx 



= -Xjr { x,y)\v-vA*f(y^mx, V )dv*dvdydx. 
We summarize the last three equalities in the following lemma. 

Lemma 4.2. Let f be a classical solution of k3.13\) subject to compactly supported initial 
conditions f$. Then, the decay of the energy functional A[f (t)] in j^.6p is governed by 

d f 

(4.7) — A[/(i)] = -A / r(x,y)\v - v*\ 2 f(y,v*)f(x,v)dv*dvdydx. 

Equipped with lemma 14.21 we can state the main result of this section. 

Theorem 4.1. Let f be the classical kinetic solution constructed in Theorem \3.1\ Then, 
the decay of its energy fluctuations around the mean bulk velocity u c , is given by 

{ C 3 e-^ , 0</?<±, 
(4-8) A[f(t)} < A[/ ] x { 

{ (i +*)-*>, a = \. 

The constants involved are = 1/(3k 3 )^(1 - A(5) > and k 5 = 2XA/^3^ > 0. 

Proof. Lemma 14. 21 implies 
d f 

—A[f(t)] = -A / r(x,y)\v - v*\ 2 f(y,v*)f(x,v)dv*dvdydx 



dt 

< -X(p(t) I \v - v*\ A f(y,v*)f(x,v)dv*dvdydx. = -2Xtp(t)M A[f(t)}. 



As before, the identity \v — v*\ 2 = \v — u c \ 2 + \ v* — u c \ 2 + 2(y — u c ) • (v* — u c ) induces the 
corresponding decomposition of the integrand on the right. Noting that 

(v — u c )f(x, v, t)dvdx = 0. 

We conclude 

|A[/(t)] < -2X<p(t)M A[f(t)] 



and Gronwall's integration yields 



(4.9) A[/(t)] < A[/oK 2A *« *(*)=/ <fi(s)ds 



o 



We distinguish between two cases. 

Case 1 [0 < (3 < |]. According to Lemma |4~T1 

cp(t) > Ak~P(i + t 2 + 1 4 )-^ > A^yPr^ \ t > 1. 
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We compute for t > 1, 

-»MS-/rtr)<lr<- ,f (t-«-l) 



which yields the first part of (|4,8p . 
Case 2 [/? = ±]. For t > 1 we have 



ex P (-2A^))<exp(-^j) <(! + *)" 



and the second part of ()4.8|) follows. □ 

5. From kinetic to hydrodynamic description of flocking 

In this section we discuss the hydrodynamic description for flocking, which is formally 
obtained by taking moments of the kinetic model (13,130 . 

(5.1a) d t f + v V x f) + AV„ • Q(f, f) = 0, x,v e R d , t > 0, 

(5.1b) Q(f,f)(x,v,t) = r(x,y)(v* - v)f(x,v,t)f(y,v*,t)dv*dy. 



We first set hydrodynamic variables: the mass p := L d fdv, the momentum, pu := L d vfdv, 
and the energy, pE 1 := 1/2 J^ d \v\ 2 fdv, which is the sum of kinetic and internal energies 
(corresponding to the first two terms in the decomposition of kinetic velocities \v\ 2 = \u\ 2 + 
\v — u\ 2 + 2(v — u) ■ u) 

1 If 

(5.2) pE = pe + -p\u\ 2 , pe:=- \v-u(x)\ 2 f(x,v,t)dv. 

For notational simplicity, we suppress time-dependence, denoting p(x) = p(x,t), u(x) = 
u(x, t) and E(x) = E(x, t) when the context is clear. 

We compute the ^-moments of f|5.1|> : multiply (|5.ip against l,v and |v| 2 /2 and integrate 
over the velocity space We end up with the system of equations, 

(5.3a) d t p + V x ■ (pu) = 0, 

(5.3b) d t (pu) + V x - (pu®u + P\ =5 (1) , 

(5.3c) d t (pE) + V x ■ (pEu + Pu + q\ = S (2) . 

Here, S^\j = 1,2, are the nonlocal source terms given by 

(,3d) 5«(M) : -A / r<x, „)(»(*) 



(5.3e) S«(M) := "A/ r (x, !/ )[ £W+E ( 9 )-„ W .„ W lp(xW 9 )* 
and P = (pij), q = (qi) denote, respectively, the stress tensor and heat flux vector, 
(5.3f) Pij ■= / (vi - Ui)(vj - uj)fdv, qi\= / (vi - Ui)\v - u\ 2 fdv. 
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Remark 5.1. The total mass of the source term <SW vanishes: exchange of variables x *-* y 
yields 

(5.4a) / S m {x , t)dx = -A f r( x, S )(»(x) - = 0. 

T/ie source term is non-positive: using 15. 2\) we find 

(5.4b) S®(t) = -\[ r(x,y)\l\u(x)-u(y)\ 2 + e(x) + e(y)}p(x)p(y)dy<0. 

JM d lz J 

We conclude that the total mass and momentum, J p(x,t)dx and J p(x,t)u(x,t)dx, are 
conserved in time. The total energy, however, 

£ (f) = / p(x,t)E(x,t)dc = - j \v\ 2 f(x,v,t)dvdx, 

JR d 2 J M 2d 

is dissipating, which is responsible for the formation of time-asymptotic flocking behavior. 
We turn to quantify this decay. We first write the total energy as the sum of total kinetic 
and potential energies, corresponding to (|5.2p . 



£(t)=£ k (t)+£ p (t), £ k (t):=lf p(x,t)\u(x,t)\ 2 dx, £ p (t) := \ [ \v-u(x,t)\ 2 fdvdx. 

2 JR d 1 JR 2d 

Lemma 5.1. The time evolution of the total, kinetic and internal energies is governed by 

d f r 1 

(5.5a) — £(t) = -A / r(x,y) -\u(x) - u{y)\ 2 + e(x) + e(y) p(x)p(y)dydx; 
at j R 2d 12 J 

(5.5b) p p(t) — — A [ r(x,y)(e(x) + e(y) - u(x) ■ - 2 / (V, . »)^; 

at J R 2d J Rd 

(5.5c) — £tM) = -- / r(x,y)( \u(x)\ 2 + \u(y)\ 2 ) p(x)p(y)dxdy + 2 / (V x -u)pedx. 
dt 2 J R 2d \ J J R d 

Proof. The equality (|5.5ap follows from integration of (|5.3cj) and invoking (|5.4bp . 
For the decay rate of the total internal energy £p (t) in (|5.5b>[) , we calculate 

d f / \ v — i/l ^ \ f \v — u\ 2 

-£p(t) = / d t r- ' ) fdvdx + / 1 - 1 d t fdvdx=:}C 1 + K2. 

dt J R 2d V 2 / J K 2d 2 

We estimate /Q separately. The first term, IC\ vanishes, for 

K\ = / {u — v) ■ dtufdvdx = I u ■ dtufdvdx — / dtu ■ (vf)dvdx 



u ■ dtupdx — I dtu ■ (pu)dx = 0. 

jR d 

For the second term, IC2, we use the kinetic model to find 

fC 2 = - [ \v - u\ 2 d t fdvdx = -- I \v - u\ 2 (v ■ V z / + XV v ■ Q(f, f))dvdx 

2 J R 2d 2 J R 2d \ / 



- f \v — u\ 2 v ■ V x fdvdx + X [ (v — u) ■ Q(f, f)dvdx =: KL21 + A/C 
>i amounts to 

K21 = -~ [ \v - u\ 2 v ■ V x fdvdx = - [ (v x \v - u\ 2 ) ■ (vf)dvdx 
2 J R 2d 2 J R 2d V / 



2 

The term /C21 amounts to 
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= — / (v — u) ■ ( (V a . • u)vf\dvdx = —2 1 (V^ • u)(pe)dx. 

jR2d V / J R d 

A lengthy calculation shows that the remaining term, IC22, equals 

IC22 = [ r(x,y)(y — u(x)) ■ (v* — v)f(x,v)f(y,v*)dv*dvdydx 
r(x, y)v • («* — v)f(x, v)f(y, v*)dv*dvdydx 
r(x,y)u(x) • (t>* — v)f{x 1 v)f(y,v*)dv*dvdydx 
r(x,y)\v — v*\ 2 f(x,v)f(y,v*)dv*dvdydx 
+7; / r(x,y)(\u{x) - u{y)f)p(x)p{y)dydx 

2 J R 2d \ J 

r(x,y)[E(x) + E{y) - -\u(x) - u(y)\ 2 ) p{x)p{y)dydx 
\ z / 



r(x, y) e(x) + e(y) + u(x) ■ u(y) p(x)p(y)dydx 
\ J 

Finally ()5.5c[) follows by subtracting (j5.5h>[) from (|5.5ap . □ 



Next, we present a fundamental estimate for the flocking behavior to the system (15.3|) . 
We set 

(5.6) r W :=/ (W^)|« + «M + 

J R 2d \2 J 

The functional T(t) can be expressed in terms of the moments Mi in (|3.9p (corresponding to 
the splitting of its integrand \\u{x) — u(y)\ 2 +e(x)+e(y) = p(x)E(x)+p(y)E(y)—u(x)-u(y)), 

T(t) = 2£(t)M - \M!\ 2 , £{t) = Mi(t) = 2£{t). 

Since Mi,i = 0, 1 are constants, this reveals that T(t) is essentially the total energy. We 
arrive at the main theorem of this section. 

Theorem 5.1. Assume (p,u,e) is a smooth solution of the system (|5.3p . (p,u,e) E C 1 ([0, T)x 
M. d ). Then we have 

r(t) < r(0)e^ 2 ^ oA *W, $(*)= /%(*)ds, <p(8):= inf r(a?(a),i/(s))cfe. 

Jo (»o,l/o) 

.Here, .A/fo is £/ie initial total mass, Mo = WpoWl 1 an d the infimum is taken over all particle 
trajectories, (x ,y Q ) \-> (x(s),y(s)). 

Proof. We use (|53al) and the relation T(t) = 2£(t)M Q - \Mi\ 2 to find 
d 



(5.7) —T(it) = -XM r{x,y)(\u{x)-u{y)\ 2 + 2e{x) + 2e(y))p(x)p(y)dydx 
at J R 2d V / 

< -XM (p(t) I (\u(x) - u(y)\ 2 + 2e(x) + 2e(y))p(x)p(y)dydx 



= -2XM <p(t)T(t). 

Gronwall's inequality then yields the desired result. □ 
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We conclude that whenever the interparticle interaction, (f(s) = inf r(x(s), y(s)) decays 
slowly enough so that its primitive, <3?(i), diverges, then flocking occurs, \u(x,t) — u(y, t)\ — > 
0, in agreement with the flocking behavior of the C-S particle model, consult remark [2T2l It is 
remarkable that the emergence of flocking is deduced here independently of the constitutive 
relation for P. In this context we observe that energy dissipation, driven by the negative 
source term in <\5.3ch vanishes as t — > oo. Indeed, theorem 15. II tells us that by (|5.4b[) . 

|5 (2) | < AT{t) -» 0. 

6. Epilogue: flocking dissipation and entropy 

We have seen that the self-propelled flocking dynamics is driven by energy dissipation. 
The dissipation mechanism reveals itself through energy decay in the particle description 
(|2.7p . in the kinetic description (|3.8bj) and equivalently, in the hydrodynamic description 
(pTTjh Observe that, 

(6.8) r(t) =M 2 M - \Mi\ 2 > 0. 

The right of (|6.8p is the usual Cauchy-Schwartz inequality 

vf(x,v)dvdx < / \v\ 2 f(x, v )dv dx x / f(x,v)dvdx. 

Thus, by theorem 15.11 <l>(i) — > oo implies time asymptotic flocking by letting T(t) — > 
which in turn, enforces an approximate Cauchy-Schwartz equality. It then follows that v 
approaches the bulk velocity, v — ► u c as t — > oo. We refer to this mechanism as flocking 
dissipation. It is intimately related to the entropy increase in the kinetic model (|3.6p . To 
this end we compute 

|/ / log(/)<fe<fo = -A/ r( x, 9 )V„lo g (/).(„-„,,//,^^ 
at J R 2d J R 4d 

r(x,y)V v f ■ (v — v*)f*dv*dvdydx 



= X r(x,y)f f*dv*dvdydx = A / r(x,y)p(x)p(y)dydx. 

This is a reversed //-theorem. Entropy increases due to the "improbable" statistical behav- 
ior of particles with increasingly highly correlated velocities, as they flock towards particle- 
independent bulk velocity. 
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